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ABSTRACT 

Millihertz quasi-periodic oscillations reported in three neutron-star low mass X-ray binaries have 
been suggested to be a mode of marginally stable nuclear burning on the neutron star surface. In 
this Letter, we show that close to the transition between the island and the banana state, 4U 1636-53 
shows mHz QPOs whose frequency systematically decreases with time until the oscillations disappear 
and a Type I X-ray burst occurs. There is a strong correlation between the QPO frequency v and 
the occurrence of X-ray bursts: when ^ 9 mHz no bursts occur, while ^ 9 mHz does allow the 
occurrence of bursts. The mHz QPO frequency constitutes the first identified observable that can be 
used to predict the occurrence of X-ray bursts. If a systematic frequency drift occurs, then a burst 
happens within a few kilo-seconds after v drops below 9 mHz. This observational result confirms 
that the mHz QPO phenomenon is intimately related with the processes that lead to a thermonuclear 
burst. 

Subject headings: X-rays: bursts. X-rays: individual (4U 1636-53, 4U 1608-52, Aql X-1) 



1. INTRODUCTION 

iRevnivtsev et all (|2001l ) discovered a new class of low 
frequency quasi-periodic oscillation (QPO) in three neu- 
tron star X-ray binary sources (Aql X-1, 4U 1608-52 
and 4U 1636-53). This new QPO has frequencies be- 
tween 7 and 9 xlO^"^ Hz, i.e. it is in the milli- Hertz 
range, and its other properties also differ from those of 
the other QPOs found i n the neutron star systems (e.g. 
energy dependeii c e, see Ivan der Klisll200i ). Although 
IRevnivtsev et al.l ()2001[ ) could not discard an interpreta- 
tion related with disk instabilities, they concluded that 
the mHz QPO is likely due to a special mode of nuclear 
burning on the neutron star surface. This interpretation 
was s trengthened by the results of lYu fc van der Klii 
(|2002t l. who showed that the kHz QPO frequency is anti- 
correlated with the luminosity variations during the mHz 
oscillation, suggesting that the inner edge of the disk 
slightly moves outward as the luminosity increases dur- 
ing each mHz cycle due to stresses generated by radia- 
tion coming from the the neutron star surface. This is 
contrary to the correlation observed between X-ray lu- 
minosity {Lx) and kHz QPO frequency, where the inner 
disk edge is thought to n iove in, as the acc retion rate 
and hence L^, increases (|van der KlisI [20061 . and refer- 
ences within). 

The properties of the mHz QPOs as observed up to 
now can be summarized as follows: (1) the fractional rms 
amplitude strongly decreases with energy, from « 2% at 
2.5 keV, down to an almost undetectable < 0.2% at w 5 
keV; (2) mHz QPOs occur only in a particular range of 
X-ray luminosity: L2-20 keV ~ 5— llx 10^^ erg/s; (3) the 
frequency of the mHz QPOs is between 7 and 9 mHz; (4) 
the mHz QPOs disappear with the occurrence of a type I 
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X-ray burst; (5) as noted above, the kHz QPO frequency 
is approximately anti-correlated with the 2 — 5 keV count 
r ate variations that c onstitute the mHz oscillation. 

iHeger et all ()2007t) suggested that the mHz QPOs 
could be explained as the consequence of marginally sta- 
ble nuclear burning on the neutron star surface. They 
found an oscillatory mode of burning, with a period Pose 
close to the geometric mean of the thermal^ and accre- 
tion'* timescales of the burning layer. For typical param- 
eters. Pose = y/Uhermai ' taeer ~ 2 miuutcs, in accordauce 
with the characteristic frequency of the mHz QPOs. The 
burning is oscillatory only close to the boundary between 
stable burning and unstable burning (in Type I X-ray 
bursts), explaining the observation that the mHz QPOs 
were seen within a narrow range of luminosities; 

Two of the three sources in which IRevnivtsev et aP 
(|2001h found the mHz QPOs are transient atoll sources 
(Aql X-1 and 4U 1608-52) whil e the third one, 4U 1636- 
53, is a persistent atoll source (jHasinger fc van der Klij 
Il989f) . The object of our current s tudy, 4U 1636-53 
has a n orbital period of « 3.8 hours (|van Paradiis et al] 
Il990f ) and a companion sta r with a mass of « 0.4 Mq 
(assuming a NS of 1.4 Mq; IGiles et al.ll2iD0l. 4U 1636- 
53 is an X-ray burst source ([Hoffman et all UWffh 
showing asymptotic burst oscillation frequencies of « 
581 Hz ([Zhang et all [19971: [Strohmaver fc Markward^ 
[2002f ). The aperiodic timing behavior of 4U 1636-53 has 
been s tudied with the EXOSAT Medium Energy instru- 
ment ([Prins fc van der Klisl[199"7l ) and with the Rossi X- 
ray Timing Explore r (RXTE, e.g. [Wiinands et al.|[1997[ : 
iDi Salvo e t"an[2003l : [Altamirano et al.|[2007[ ). 

4U 1636-53 is a reference source for studying nuclear 
burning on the surface of a neutron star since it shows 
the full range of burst behavior: single and multi-peaked 

^ The thermal timescale is defined as tthermal = (^pT/t where Cp, 
T and e are the heat capacity at constant pressure , the temperature 
and the nuclear energy generation rate, respectively 

* The accretion timescale is defined as taccr = y/rh where y and 
rh are the column depth of the burning layer and local accretion 
rate, respectively. 
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Fig. 1. — Color-color diagram as lAltamirano et all 1)20071 ). 
Each data point represents the average of an observation (Ri 2 
to 30 ksec). The ellipse marks the region in which mHz QPOs 
with decreasing frequency were found. The labels A, B and C 
correspond to those in Figure |3] 
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Fig. 2. — Dynamical power spectrum smoothed with a 750 sec- 
onds sliding window with steps of 200 seconds, showing the mHz 
QPOs during the last 12 ksec before the X-ray burst occurs. This 
sequence corresponds to case B in Figures |3] & [U — Obsid: 60032- 
05-02-00. The three black vertical lines correspond to the times of 
occurrence of the X-ray bursts. For clarity, we plot only powers 
above 10 which correspond to > 2a (single trial per 750 seconds 
window but normalized to the number of possible frequencies in 
the range 0-0.5 Hz). 

Type I X-ray bursts, superbursts, burst oscillations, pho- 
tospheric radius expansion, regular and irregular burst 



sequences (e.g. lGallowavet"aIll2006D and mHz QPOs. 
As such, it is an ideal source to understand the relation 
between these different observational manifestations of 

nuclear bur ning. 

Recently. IShih et all (|2005l ) reported that 4U 1636-53 
has shown a significant decreas e in its persistent dur 



ing the years 2000 and 2001. lAltamirano et al.l |2007f ) 
show that during the low period, 4U 1636-53 is ob- 
served in its (hard) island states. This provides an op- 
portunity to study the mHz QPOs in harder and lower 
luminosity states than was possible up to now. 

2. DATA ANALYSIS & RESULTS 

We used data from the RXTE Proportional Counter 
Array and the High Energy X-ray Timing Experiment 
(PCA a nd HEXTE, respectively; for instrument informa- 
tion see iJahoda et ani2006l : iGruber et al.lll996l ). Up to 
June, 2006, there were 338 pubhc pointed observations. 
An observation covers 1 to 5 consecutive 90-min satellite 
orbits. Usually, an orbit contains between 1 and 5 ksec 
of useful data separated by 1-4 ksec data gaps; on rare 
occasions the visibility windows were such that RXTE 
continuously observed the source for up to 27 ksec. In 
total there were 649 gap-free data segments of length 0.3 
to 27 ksec. 

We produced energy spectra for each observation using 
Standard data modes and fitted them in the 2 — 25 keV 
and 20 - 150 keV bands for PCA and HEXTE, re- 
spectively. The interstella r absorption T Vg was fixed 
at 3.75 X 1021 cm-2 (see iSchuld [19991 : iFiocchi et all 
[2006h . We used 1 -sec resolution event mode PCA light 
curves in the « 2 — 5 keV range (where the mHz QPOs 
are strongest) and searched for periodicities in each of 
the 649 se gments separately using Lomb-Scargle per i- 
odograms (|Lomblll976l: IScargldfl98a iPress et al.l[l99l . 
Segments in which one or more Type I X-ray bursts were 
detected were searched for periodicities before, in be- 
tween, and after the bursts. We find that the oscillations 
in the w 2 — 5 keV range ar e evident from the light curves 
(see for example figure 1 in lRevnivtsev et al.|[200ll ). The 
significanc e as estimated fr om our Lomb-Scargle peri- 
odograms ijPress et al.lll992l ) confirm that the oscillations 
are all above the 3a level. We estimated the uncertain- 
ties in the measured frequencies by fitting a sinusoid to 
1000 sec data segments to minimize frequency-drift ef- 
fects. The typical errors on the frequency are of the 
order of 2 - 6 x 10"^ Hz (or 2 - 6 x 10"^ mHz). 

We detected mHz QPOs in 124 of the 649 segments. 
Most occur in segments with less than 4 ksec of useful 
data and sometime s the QPO s cover only part of a seg- 
ment. Revnivtsev e t al.l (|2001h reported the characteris- 
tics of the mHz QPOs between March 1996 and February 
1999. Using the X-ray colors average d per observation as 
reported bv lAltamirano et "al] (|2007( ). we find that their 
data sample the region at hard colors < 0.7 and soft col- 
ors > 1 (see Figure IH , which r epresents the so called 
banana state (Ivan derKlisl[200 6!). Some of the later ob- 
servations also sample the banana state. We re-analyzed 
all the data in this region of the color-color diagram and 
fo und results which are co nsistent with those reported 
bv iRevnivtsev et al.l ()2001[ ): the frequency of the QPOs 
varies randomly between 6 and 9 mHz. 

In the harder state close to the transition between the 
island and banana state and marked with the ellipse 
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Fig. 3.— 2-60 keV PCA light curve of observations 60032-05-01,02,03 and 04. A, B and C indicate intervals in which the mHz QPOs 
are detected. In all three cases, the frequency of the niHz QPO decreases with time until the appearance of an X-ray burst. 



drawn in Figure [U we found 22 segments with signifi- 
cant mHz QPOs; in these observations, the (2-150 keV) 
luminosity was 6 — 10 x [d / {Qkpc))^ erg s~^, while 
for the other cases, corresponding to the banana state, it 
was higher (10 - 35 x W%dl{<okpc)f erg s.-^). 

Among the 22 segments, we distinguish two groups 
based on segment length: the first consisting of 4 seg- 
ments each with more than 14 ksec of uninterrupted 
data, and the second consisting of 18 segments each cor- 
responding to one orbit with less than « 4 ksec of useful 
data. For all four segments in the first group, we mea- 
sure a systematic decrease of frequency from between 
10.7 and 12.5 mHz down to less than 9 mHz over a time 
interval of 8 to 12 ksec, after which an X-ray burst oc- 
curs and the QPOs disappear (the QPOs become much 
less than Scr significant). Figure [2] shows a representative 
dynamical power spectrum corresponding to one of these 
segments (interval B in Figure [3]). The QPO is present 
ft! 12 ksec seconds before the burst and its frequency sys- 
tematically decreases with time from « 10.7 mHz down 
to « 7.6 mHz. Then the X-ray burst occurs and the QPO 
disappears. In the second group, 16 of the 18 segments of 
< 4 ksec show the mHz QPO frequency to decrease either 
within a segment, or between 2 or 3 consecutive orbits 
(with 2-4 ksec data gaps in between) at rates consistent 
with those seen in the 4 long segments. The two remain- 
ing segments are too short and isolated to constrain the 
frequency drift well. 

To illustrate the interplay between this very system- 
atic behavior of the mHz QPOs and our data structure, 
in Figure[3]we show a representative lightcurve. A, B and 
C mark three intervals in which the mHz QPOs were de- 
tected and that each terminate with an X-ray burst. As 
can be seen, we have data in which mHz QPOs are de- 
tected and followed through consecutive segments (case 
A), and data in which the oscillations are detected and 
disappear within one segment (cases B & C). Further- 
more, we have data in which the oscillations are present 
from the start of the observation (case B) as well as data 
in which the mHz QPOs appear during an observation 
(case A & C). 

Among our 22 segments, the frequency of the oscil- 
lations varies in the range 7 — 14.3 mHz with directly 
observed onset frequencies between 10.7 and 14.3 mHz. 
Interpolating through gaps, the QPOs last for 7.5 to 16 



ksec. Over such intervals, the frequency is always consis- 
tent with decreasing at average rates from 0.07 to 0.15 
mHz ksec~^, and the frequency had always dropped to 
< 9 mHz just before an X-ray burst (as estimated from 
the last 750 seconds before the burst). Interestingly, this 
last result applies to all cases in which we detect the mHz 
QPOs before an X-ray burst, including the cases that oc- 
cur in the banana state: it seems that independent of the 
spectral state of the source, no X-ray burst will occur if 
the mHz QPOs are present at a frequency higher than 
» 9 mHz (bursts do occur in both states that are not 
preceded by detectable mHz QPOs). 

No relation between the 2-60 keV count rate and fre- 
quency was found: in two of the four long segments the 
count rate decreased about 10% during the time the mHz 
QPOs were present, while in the other two cases the 
count rate increased by approximately the same amount. 
No clear relation was found between frequency range cov- 
ered and duration of the oscillation; perhaps this is re- 
lated to the fact that, as shown in Figure [21 the fre- 
quency does not decrease smoothly but has short periods 
in which it is consistent with being constant. 

When 4U 1636-53 is observed close to its island- 
banana state transition, the mHz QPOs disappear only 
when an X-ray burst occurs. However, this is not the case 
for the banana state, in which we found also observations 
in which the oscillation disappear below detectable levels 
without the occurrence of an X-ray burst. The interval 
of time required to again detect the oscillations after a 
burst occurred is variable. The two extreme cases are (i) 
observation 60032-01-06-000 where no mHz QPOs were 
detected during the « 15 ksec of uninterrupted data fol- 
lowing an X-ray burst and (ii) observation 40028-01-06- 
00, where mHz QPOs are detected again k, 6000 seconds 
after an X-ray burst occurred. We note that in the first 
case, the source was close to the transition between island 
and banana state while in the second case, the source was 
in the banana state. Nevertheless, no clear relation be- 
tween this waiting time and the source state (island or 
banana state) was found. As bursts may be missed due 
to data gaps, a time interval of « 1000 seconds between 
a (missed) X-ray burst and onset of the QPOs cannot in 
some cases be excluded (e.g. case A in Figure [3]). 

3. DISCUSSION 
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We have shown that close to the transition between 
the island and the banana state 4U 1636-53 exhibits 
mHz QPOs whose frequency systematically decreases 
with time until the oscillations disappear with the occur- 
rence of a Type I X-ray burst. The mHz QPO frequency 

V constitutes the first identified observable that can be 
used to predict the occurrence of X-ray bursts: when 

> 9 mHz no bursts occur, while v mHz does allow 
the occurrence of bursts. If a systematic frequency drift 
occurs, then a burst happens within a few kilo-seconds 
after v drops below 9 mHz. This observational result 
confirms that the mHz QPO phenomenon is intimately 
related with the processes that lead to a thermonuclear 
burst. 

The fact that the observation of a systematic frequency 
decrease with time implies the occurrence of a future X- 
ray burst, strongly suggests that the frequency of the 
mHz QPOs is related to the burning processes on the 
neutron star surface. One possibility is that the fre- 
quency of the QPO is somehow a measurement of the ac- 
cumulation of fresh fuel on the neutron star surface which 
will be available for a future thermonuclear burst. To our 
knowledge, there has been only one attemp t to theoret- 
ically explain the mHz QPOs phenomena (jHeger et al.l 
|2007() . In this model the frequency of the QPO depends, 
among others, on the amount of available fresh fuel, on 
the local accretion rate and the composition of the ma- 
terial. It is beyond the scope of this Lette r to perform 
numer ical simulations as those reported bv iHeger et al.l 
(|2007| ). In the rest of this discussion we briefly compare 
these authors' model predictions with our observations 
and propose some more complex scenarios. 

Analytical and numerical results bas e d on the sim- 
plified one-zone model of iPaczvnskil (jl983f ) in the 
IHeger et al.l (l200l marginally stable burning model (see 
Section [1]) predict that (i) close to the boundary be- 
tween stable and unstable burning, the NS surface will 
show temperature fluctuations with constant frequency 

V if the local accretion rate m remains constant; (ii) 
this marginally stable burning regime will occur at m 
near Eddington, hence accretion must be confined to 
a surface area Sa that is much smaller that the total 
ar ea of the NS; (iii) v correlates with m (see figure 4 
in IHeger et al.ll2007| ) and (iv) thermonuclear bursts and 
mHz QPOs should not be observed at the same luminos- 
ity and therefore presumably at the same m. 

In this paper, we show that for constant luminosity 
the QPO frequency can systematically decrease in time 
and that instantaneously measured frequencies can be 
the same for different luminosities. We also show that 
mHz QPOs and thermonuclear bursts do in fact occur at 
the same luminosity and that both phenomena are clearly 
related. This means that we are dealing with a more 
compl ex scenario than that introduced by IHeger et al.l 
(1200% . 

The amount of time between the preceding X-ray burst 
and the onset of niHz QPOs is variable (> 6ksec) and 
apparently independent of source state. If the system is 
locally accreting at m ~ m Edd and if none of the accret- 
ing fuel is burnt, only k, IQQQ seconds are required to 
accrete a fuel layer of column depth y / capable of under- 
going marginally stable burning (i// k, 10^ g cm~^ and 
m w 8 • 10** g cm"^ s~^ - see e.g. IHeger et al.ll2007f) . One 



possible explanation for the observed longer intervals be- 
tween burst and onset of oscillations, is that a large frac- 
tion of the accreted fuel is burnt as it is accreted on 
the neutron star surface. Of course the burning fraction 
could vary in time, and this estimate is assuming that all 
the fuel was bur nt during the l ast X-ray burst, which is 
not always true (|Bildstei]||1998( ). Interestingly, if this in- 
terpretation is correct and low partial burning fractions 
can occur, under certain conditions the niHz QPOs could 
appear in much less than a 1000 seconds after an X-ray 
burst. 

The fact that the amount of time between the preced- 
ing X-ray burst and the onset of mHz QPOs is variable 
may be also an indication that not all the accreted fuel 
is burnt nor available to participate in the marginally 
stable burning. For example, accretion could occur onto 
an equatorial region o ccupying less than 10% of the sur- 
face area of the star ()Heger et al.ll2007l) . A possibihty 
is that part of the fresh fuel burns stably at a rate B(t) 
per unit area while the other part leaks away from this 
region at a rate R(t). While the material accumulated 
at a rate i?(t) would serve as fuel for a thermonuclear 
burst, marginally stably burning of the matter on the 
equatorial belt is (in principle) still possible. Although 
such scenario cannot explain the frequency drifts we ob- 
serve, it can explain why mHz QPO and X-ray bursts do 
occur at the same m. If mHz QPOs can only occur at a 
certain local accretion rate m (~ rriEdd), a small change 
in effective local accretion rate will lead to an absence of 
mHz QPOs. This might explain why the mHz QPOs are 
not always present between X-ray bursts. 

Another pos s ibility (which is not taken into account 
in IHeger et al.l (|2007l )'s model) is that there is a signif- 
icant heat flux from deeper in the star that heats the 
region undergoing marginally stable burning. For exam- 
ple, changes in heat flux due to energy that is first con- 
ducted into deeper layers during an X-ray burst and then 
slowly outwards towards the surface might be possible. 
Such a change in the heat flux could affect the conditions 
of the burning layer (e.g. temperature or burning rate 
B{t)) and therefore could affect the characteristics of the 
burning processes on the neutron star surface. 

Other aspects of the observations offer further chal- 
lenges for theoretical models that explain burning pro- 
cesses on the neutron star surface as well those which 
explain atoll sources states. In particular, (i) why the 
systematic frequency drifts are observed close to the tran- 
sition between the island and the banana state while 
the frequencies are approximately constant in the ba- 
nana state. This may be another indication that the 
disk geometry of th e system is changiiig du ring the state 
transition (see e.g. iGierlihski fc Doiiil200 2') : (ii) why in 
the transition between island and banana state the oscil- 
lations disappear only when an X-ray burst occurs, while 
in the banana state they can also disappear without an 
X-ray burst (see Section [2]). Clearly, further theoretical 
work is needed. More observational work on the interac- 
tions between mHz QPOs and X-ray bursts is in progress 
and will provide further clues for theoretical models. 
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